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Abstract. The Virgo gravitational wave detector is an interferometer (ITF) with 
3km arms located in Pisa, Italy. From July to October 2010, Virgo performed its third 
science run (VSR3) in coincidence with the LIGO detectors. Despite several techniques 
adopted to isolate the interferometer from the environment, seismic noise remains an 
important issue for Virgo. Vibrations produced by the detector infrastructure (such 
as air conditioning units, water chillers/heaters, pumps) are found to affect Virgo's 
sensitivity, with the main coupling mechanisms being through beam jitter and scattered 
light processes. The Advanced Virgo (AdV) design seeks to reduce ITF couplings to 
environmental noise by having most vibration-sensitive components suspended and in- 
vacuum, as well as muffle and relocate loud machines. During the months of June and 
July 2010, a Guralp-3TD seismometer was stationed at various locations around the 
Virgo site hosting major infrastructure machines. Seismic data were examined using 
spectral and coherence analysis with seismic probes close to the detector. The primary 
aim of this study was to identify noisy machines which seismically affect the ITF 
environment and thus require mitigation attention. Analyzed machines are located at 
various distances from the experimental halls, ranging from 10m to 100m. An attempt 
is made to measure the attenuation of emitted noise at the ITF and correlate it to the 
distance from the source and to seismic attenuation models in soil. 
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1. Introduction 

The general theory of relativity predicts that all accelerating objects with non-symmetric 
mass distributions produce gravitational waves (GW) i". LIGO (the Laser Interferometer 
Gravitational- Wave Observatory) [1] and Virgo [2] experiments seek to directly detect 
GWs and use them to study astrophysical sources. 

Virgo, located in Cascina, Italy, consists of a laser Michelson interferometer with 3km 
long, Fabry Perot resonant optical cavities in its arms. With respect to other similar 
detectors, Virgo has enhanced sensitivity between 10 and 100Hz due to the seismic 
isolation performance of the "super-attenuators" to which the test masses are suspended 
[3]. Beginning in 2011, the detector expects to undergo upgrades, known as Advanced 
Virgo (AdV), to improve its sensitivity by one order of magnitude [4]. 

Seismic noise places a limit on Virgo's detection sensitivity. In particular, the main 
source of seismic noise in the 10 to 100Hz band at the Virgo site is the detector's 
infrastructure machines (please see the left of Fig. 1 for a comparison with a low 
noise model). Although Virgo's mirrors are well-isolated from local seismic activity 
by suspension systems made of multi-stage pendulums, seismic noise remains a concern. 
The residual ground motion causes "diffused light" [5]. Because of unavoidable 
imperfections in the detector's optical components, some tiny fraction of light can exit 
the main optical path and strike a surface that is connected to the ground, and thus be 
subject to the local seismic field. When this light scatters off of objects connected to 
the ground (such as optical components on external tables used for detector controls), 
it is often diffused over a wide solid angle. A fraction of the light can re-enter the main 
beam path, but with a noisy phase modulated by the seismically excited scattering 
object. This phenomenon functions as additional noise in the GW channel and limits 
its sensitivity. 

Because seismic events couple in this way into the GW channel, it is necessary to 
understand and attenuate the noise from the local seismic environment. Many of the 
strong lines seen during VSR2 can be attributed to seismic sources (please see the right 
of Fig. 1). To assist in the identification effort, the Virgo detector is supplemented with 
several types of environmental sensors, including seismometers and accelerometers, that 
monitor the local environment [6]. These channels are used to detect environmental 
disturbances that can couple to the GW channel and are placed in sensitive areas of the 
interferometer. To reduce the influence of anthropogenic noise, during the Advanced 
Virgo upgrade, machines that are identified as seismically and acoustically affecting 
the interferometer will be replaced, moved, or isolated [1]. The current proposal is to 
move all chillers, water pumps, air compressors, and air conditioners from their current 
locations to other areas further from the interferometer. By placing these machines 
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on their own isolated platforms and by adopting techniques to reduce noise emission 
and propagation, Virgo hopes to significantly reduce their effects on the detector. 
To prioritize the actions required to reduce noise at Virgo, a careful documentation 
of the present local seismic environment is necessary, as well as identification among 
the infrastructure machines of relevant sources of seismic noise reaching the detector's 
sensitive components. The work described hereafter is part of this effort. 




Frequency [Hz] '" frequent [Hz] 



Figure 1. Left: Variation in the root spectral density of the Guralp seismometer at 
the Virgo site during a day in June of 2010. The two solid curves correspond to the 
Peterson low- and high- noise models [7]. Right: Plot of the VSR2 sensitivity as well 
as the coherence projection of seismic probes on the sensitivity. 



During the months of June and July 2010, a Guralp-3TD seismometer was stationed 
to monitor floor vibrations at various technical areas around the Virgo site which host 
several infrastructure machines (i.e. water chillers, heaters, pumps). The power spectra 
of the measured seismic data were examined, and several single sources were identified 
using correlation of the seismic data with synchronous records of the infrastructure 
machines working status (i.e. temperature monitors, voltage sniffers). Then, seismic 
data were correlated with the synchronous record of permanent seismic probes close 
to the detector. Studying the calculated coherence and power spectra, we selected 
"relevant" sources whose noise emissions reach the detector and measure the noise 
attenuation between the sources and detector. We also measured approximate noise 
reduction factors as a function of the distance between the machine and the detector. 

Section 2 discusses the methods used when taking and analyzing data. Section 3 
describes examples of noise emissions from the sources found in this analysis at the 
various locations, while Section 4 provides an estimate for the attenuation of the noise 
spectral peaks as a function of distance. Section 5 presents our conclusions. 
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2.1. Measurement Setup 

We acquired several sets of data with a Guralp 3TD tri-axial seismometer [8] at 
six locations around the site. The chosen sampling rate was 250 Hz. The Guralp 
seismometer is equipped with a GPS antenna receiver to synchronize its data to Virgo 
seismic channels. Maps for the site and the locations of the measurements can be seen 
in Fig. 2. These locations were chosen for their proximity to machinery known to create 
large seismic noise. These locations are detailed in Table 1. In each location, the Guralp 
seismometer was installed on the same concrete platform where the machinery sat and if 
possible, at a distance of at least a few meters from the machinery to prevent saturation. 
For each location, the Power Spectral Density (PSD) of the Guralp seismometer, referred 
to as the "test probe," is computed (see the left of Fig. 1). Each measurement lasts 
approximately 24 hours in order to determine the hourly and daily cycle operation of 
the machines. This measurement is compared with a tri-axial seismometer, referred to 
as the "reference probe," permanently stationed in the closest nearby experimental area 
containing sensitive interferometer components [6]. In the Mode Cleaner, West End, 
and North End Buildings, the reference probe was an Episensor FBA ES-T, while in 
the Central Building, the probe was a Guralp CNG-T40. If the source is not obvious, 
suspected machines are analyzed with a piezoelectric accelerometer (PCB) [9] placed in 
direct contact with the machine and read out with a spectrum analyzer. The machines 
characteristic frequencies are then compared to those seen in the coherence. 




Figure 2. Left: Aerial View of the Virgo Interferometer. Right: Aerial view of Virgo's 
Central Area. Measurement sites are marked with letters (see Table 1). 
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Section Letter 



Test Probe Location 



Reference Probe Location 



Distance Between Seismometers (m) 



A 



Technical Building One: Ground Floor 

Building One: First Floor 
Outside Building One: Ground Floor 
Outside Mode Cleaner Building: Ground Floor 
West/North End Building: Ground Floor 



Central Building: Ground Floor 
Central Building: Ground Floor 
Mode Cleaner: Ground Floor 



Ml 



B 



105 



C 



30 
11) 



L) 



Mode Cleaner: Ground Floor 
West /North End Building: External Optical Bench 



E/F 



211 



Table 1. Locations of seismometer placement. 



2.2. Analysis Methods 

We considered relevant the frequency lines (persistent spectral peaks) that appear as 
significantly coherent (coherence rising out of the background) between the test and 
reference probes, as the source of those lines is very likely to be the same.. Coherence 
itself does not allow one to distinguish between the seismic signals generated close to the 
test probe (i.e. machinery under test) from those produced close to the reference probe 
(i.e. cooling fans of electronic units). In order to determine when the former is the case, 
the ratio of the test and reference seismometers PSDs is examined, and coherent lines 
that are stronger near the reference probe are selected. These techniques work well for 
machines producing continuous lines, such as water pumps and cooling fans. 

Periodic lines, on the other hand, are more difficult to identify. Because the lines come 
and go, they tend to be washed out in PSD averaging. For this reason, frequency-time 
(FT) plots of the PSDs are produced and examined by eye for periodic lines. In order 
to identify the source of these lines, we compute the root-mean-square (RMS) around 
that line, which is defined as: 



where Xi is the ith component of the frequency band in question and 5i is the width 
of the frequency bins in the spectrum. The sum goes from i — fx to i — /2, where /i 
and f2 are the minimum and maximum frequencies respectively in the frequency band. 
This RMS value is then correlated to the time series of various Infrastructure Machine 
Monitoring System (IMMS) signals, including temperature and pressure probes. 

3. Characterization Examples 




(1) 



In this section, the most significant examples utilizing methods discussed above are 
presented. In Table 2, information about important, identified noise lines are given. 
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Location 


Frequency (Hz) 


Periodicity (Period) 


Source 


PSD Ratio 


Technical Building One 


24.2, 48.4,78.6 


Periodic (42 Minutes) 


Cold Water Chiller 1 


16 


Building One 


19.3 


Continuous 


Computer Fans 


13 


Outside Building One 


48.9 


Continuous (Daytime) 


Bl Chiller 


2 


Mode Cleaner 


48.9, 97.8 


Periodic (21.5 Minutes) 


Mode Cleaner Chiller 


220 


West End Building 


47.1 


Continuous 


Warm Water Pump 


1.1 


West End Building 


48.8 


Continuous 


Cold Water Pump 


6.1 


North End Building 


22.8 


Continuous 


W r ater Pump 


4.2 



Table 2. Important noise lines identified in the study. The table provides the line's 
frequency, periodicity, source, and the ratio of the PSDs of the test and reference 
probes. 



3.1. Characterization of Noise from Technical Building One 

A data set was taken in Technical Building One (TBI) and compared with data from 
a seismometer on the floor of the Central Building (CB), where the two arms of the 
interferometer converge. These two probes are approximately 80 meters apart. The 
PSDs of the test and reference probes, as well as the coherence between them, can be 
seen in Fig. 3. A periodic line around 24.2Hz (with harmonics at 48.4Hz and 72.6Hz) 
is seen, which is showed in detail on the left plot of Fig. 4. 

PSD Comparison Plot Coherence Between Seismometers in B1 and CB 




Frequency [Hz] Frequency [Hz] 

Figure 3. Left: PSDs of the seismometers in TBI and the CB. The peaks topped by 
red stars correspond to the 24.2Hz periodic line while the blue stars correspond to the 
48Hz continuous line. Right: Coherence between seismometers in TBI and the CB. 
The peak topped by the red star corresponds to the 24.2Hz periodic line while the blue 
star corresponds to the 47.9Hz continuous line. 

The coherence for this line between the test and reference probes was about 0.6, 
indicating a relevant correlation between the two seismometers at that frequency. In 
order to find the source of the periodic line, it was necessary to determine in which 
building the line was louder. The ratio of the PSDs (Seismometer in TBI /Seismometer 
in CB) was approximately 16, indicating the source of the line to be in TBI. 

The time series of several machinery monitors in TBI were compared with the RMS 
of the 24.2Hz signal in the Guralp seismometer, one of which may be seen on the right 
of Fig. 4. From this figure, the RMS is clearly correlated with the time series of a 
temperature monitor of the first water chiller. We see that when the temperature of the 
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water reaches a high point, the cold water chiller switches on, causing the temperature 
of the water to decrease. When the temperature reaches a low point, the water chiller 
switches off and the whole process starts again. The chiller is located on the roof of 
TBI. Even though the chiller is equipped with insulating springs, the vibrations can 
probably travel through the rigid water pipes between TBI and the CB or in the water 
itself. 

Time-Frequency PloHor the Guralp Seismometer Centered Around 24Hz V1JMMS_TETB1 1JNLCWCHIL1 vs. Guralp 24 Hz Comparison 

Guralp 24 Hz Line 

V1JMMS_TETB11JNLCWCHIL1 . a , n , 

5 10 15 '0 12345678 

Time Time (Hours) 

Figure 4. Left: FT-plot of the test probe in TBI zoomed in on the 24 Hz region. 
Right: A plot of both the RMS in the 24Hz band of the Guralp seismometer as well 
as the time series of a temperature monitor of the first (of two) cold water chillers. 

3.2. Characterization of Noise from the Cold Water Chiller Outside the Mode Cleaner 

Another data set was taken outside the Mode Cleaner Building (MCB), on the same 
concrete platform as the MC's cold water chiller, and compared with the seismometer 
on the MCB floor. The MCB contains the end mirror of the Input Mode Cleaner, an 
optical cavity that filters jitter and power noise as well as higher order modes from the 
beam. The distance between these probes is approximately 10 meters. 

In viewing the PSD plots, a number of significant lines are noticed. This includes 
continuous lines at 24.2Hz and 38.8Hz and a strong periodic line at 48.9Hz, which has a 
harmonic at 97.8Hz. To study this line, the RMS around the 48.9Hz band was computed 
and compared to probes inside the MCB. A water temperature probe follows this line 
well, as can be seen in Fig. 5, indicating its source is the MC water chiller. This line is 
important since coherence was noted at its frequency between the GW channel of VSR2 
run data and the seismometer in the MCB (see right of Fig. 1). Thus, it seems that 
the MCB water chiller affects the Virgo sensitivity and for this reason, requires prompt 
attention for noise reduction. 

3. 3. Characterization of Noise in the West End Building 

Another data set was taken inside the West End Building (WEB) near the machinery 
and compared with an Episensor seismometer located on the external optical bench. 
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VI JMMS_TEMC51 JNLCW vs. Guralp 49 Hz Comparison 



1 o 



— Guralp 49 Hz Line 

— V1JMMS_TEMC51JNLCW 




1.5 2 2.5 
Time (Hours) 



Figure 5. Comparison between the RMS in the 48.9Hz band of the Guralp 
seismometer and the time series of a temperature monitor of the water chiller. 



The distance between these seismic probes is approximately 20 meters. 

Analyzing the PSD plots, we noticed several significant lines, including continuous 
lines at 24.7Hz, 47.1Hz, and 48.8Hz. The 47.1Hz and 48.8Hz are associated with the 
"warm" and "cold" water pumps respectively as verified with the PCB probe in Fig. 6. 
As the 47.1 Hz signal reached the optical bench with almost no attenuation, there is 
likely a preferred path, for example through water pipes, which must be investigated 
further. 



PSD Comparison Plol 



-Warm Water Pump 
-Cold Water Pump 




30 40 50 
Frequency [Hz] 



Figure 6. This plot shows the PSD for both the warm and cold water pumps in the 
WEB. The data were taken with the PCB probe in direct contact with the machines 
under study. The warm water pump's PSD has a strong line around 47.1Hz, which is 
shown by the red star on the left. The cold water pump's PSD exhibits a strong line 
around 48.8Hz, which is shown by the blue star on the right. 



4. Attenuation Measurement 



Using the results above, we seek to measure how much noise attenuation occurs when 
a machine is on its own foundation and a certain distance away from an experimental 
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area. The PSD ratio of test and reference probes is used as a tentative estimate of the 
attenuation of the vibration, the result for which can be seen in Fig. 7. Measurements are 
compared to a dissipation law accounting for geometrical spreading of surface cylindrical 
waves and for energy dissipation in soil [10]. 

As can be seen in the plot, distances of about 30 meters or less have an attenuation 
of approximately a factor of 3, while distances of about 80 meters and greater have 
attenuation values greater than 10. Indeed, the simple noise attenuation model does 
not apply in most of the examined cases. One reason is that the point-like source 
approximation is not valid for short propagation distances (d < A, with our seismic waves 
length ~ 10m) or for extended sources (L > d, with the machine platform size L ~ 10m). 
In addition, the soil medium is often not uniform, and waves' reflection and refraction 
can occur at discontinuities. Wave amplification due to mechanical resonances of the 
two platforms (technical area and experimental area) which are excited by the seismic 
noise is also possible. Another important reason is that in most cases, propagation does 
not occur uniquely through soil. Less dissipative paths can exist, such as through water 
pipes or pressure waves inside the water itself. 
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Figure 7. A plot of the average of the PSD ratio of the coherent lines as a function 
of distance between the probes. We selected locations separated by at least a few 
wavelengths of soil. The attenuation factors data at 105m are lower limits only. 

5. Conclusion 

Noise from several Virgo infrastructure devices, such as water chillers, heaters, and 
pumps, seismically affects sensitive parts of the interferometer. Although an impact on 
the present interferometer is not evident in Virgo's sensitivity (with the exception of the 
MC chiller), the noise reaching the experimental area is, in several cases, considerably 
above the background, and requires a reduction for AdV. 
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Measured attenuation of seismic signal from a distant source seems not compatible 
with the much stronger attenuation expected from soil dissipation. It is suspected that 
water pipes might function as seismic shortcuts and thus need more comprehensive 
investigation and mitigation attention. Noise transmitted through soil can be more 
efficiently reduced using seismic isolation systems (i.e. springs). 

Further beneficial studies might include setting up a noise source with known power and 
characteristic frequencies at various distances from a sensitive experimental area. This 
would allow a comprehensive study of the noise attenuation as a function of distance, 
something we were unable to do in this study. 
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